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The vulnerability to real-life networks against small initial attacks has been one of outstanding challenges in 
the study of interrelated networks. We study cascading failures in two interrelated networks S and B 
composed from dependency chains and connectivity links respectively. This work proposes a realistic model 
for cascading failures based on the redistribution of traffic flow. We study the Barabasi- Albert networks 
(BA) and Erdos-Renyi graphs (ER) with such structure, and found that the efficiency sharply decreases with 
increasing percentages of the dependency nodes for removing a node randomly. Furthermore, we study the 
robustness of interrelated traffic networks, especially the subway and bus network in Beijing. By analyzing 
different attacking strategies, we uncover that the efficiency of the city traffic system has a non- equilibrium 
phase transition at low capacity of the networks. This explains why the pressure of the traffic overload is 
relaxed by singly increasing the number of small buses during rush hours. We also found that the increment 
of some buses may release traffic jam caused by removing a node of the bus network randomly if the damage 
is limited. However, the efficiencies to transfer people flow will sharper increase when the capacity of the 
subway network > ao. 



Many real- world networks interact with and depend upon other networks \ In the past few years, inter- 
connected networks^"^ that featured in the interaction among the constituents of a variety of complex 
systems referred to biological, technological and social systems. The researchers shifted their attention to 
the multi-layered^"^ character of real-world systems that were coupled and especially interdependent with many 
networks. Real systems are studied to be composed by rather small sub-groups in which the nodes belong to a 
chain strongly depended on each other^. In our daily life, we experience social behaviors by distinguishing various 
roles within each single layer we belong to. A network was proposed for two interdependent networks and two 
types of link containing: connectivity and dependency links^°'^\ Take the city traffic system for example, the 
subway network and the bus network constitute the daily traffic network, all nodes in different subway lines 
constitute the subway network, and the nodes of the same subway line depend on each other. Recent studies have 
shown that the topology of each interdependent layer can affect dramatically the properties of percolation^^ and 
cascade failures^^ in dynamic complex systems. Such systems were found to be more vulnerable compared to 
classical networks with only connectivity links. Their percolation transition is usually a first-order one compared 
to the second-order transition found in classical networks^^'^^. City traffic networks as one of such systems, which 
couple bus and subway networks, has caused much attention^^'^^. Cascading failures are common in most of the 
coupled communication and transportation networks. The failure of initial removal nodes in one network will 
cause further failures in other networks and vice versa^^ resulting in a cascade of failures. Buldyrev et al.^° studied 
the robustness of two interacting networks, and found that interdependent networks became significantly more 
vulnerable compared to their non-interacting counterparts. Cao et al.^^'^° analyzed the robustness of a network of 
networks, and gave a general analytical framework for studying percolation of n interdependent networks. It was 
found by Schneider et al.^^ that scale-free networks with an onion-like structure are very robust against targeted 
high degree attacks. An optimal way of interconnection was found against cascading failures^^'^^ Recently, these 
works of Refs. [9,11,24,25] demonstrated that a crossover in phase transitions of cascading failures can arise in 
interdependent networks. 

The unprecedented pressure of traffic congestion forced many researchers to study the robustness of city traffic 
networks. Understanding the robustness of city traffic networks has become a major challenge for solving the 
traffic congestion problem. In previous studies, the failure nodes only cause the failure of nodes that are connected 
to them directly in one or more networks. The interrelated traffic systems contain dependency subway chains and 
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Figure 1 | Cascading failures in the dependent network (5) and the connected network (B) with the initial removal of a single node chosen at random in 
the connected network B. The colorful solid spheres are different sites. Nodes 1, 2 and 3 linked with broken lines in network S are in the dependent chain 
1. Nodes 4 and 5 being connected with solid lines in network S belong to dependent chain 2. The nodes with the same number in these two networks are 
the nodes that have both dependent and connected edges. Suppose that node 3 in network B is attacked at t = 0, so nodes 3 in network S also could not 
work at t = 1. Then nodes 1, 2 and 3 of network S dependent linking with network S could not work in the next time step t = 2. This may cause the 
congestion of some relevant nodes in the connected network and an iterative cascading of overload between both networks. 



connectivity nodes of buses. A subway network S is composed of 
nodes, and the subway lines often intersect at the same nodes. The 
subway network has a special property, i.e. if a node of a dependent 
subway line (or lines) failed, then it will lead to the failure of its entire 
dependent nodes of the same line (or lines). That in turn leads to a 
new failure among connecting nodes of the bus network B. Once the 
cascade process is triggered, it will only stop if nodes that fail in one 
step do not cause additional failure in the next step. The existing 
models cannot explain this appropriately. Therefore, to understand 
the intricate cascade failures of real complex traffic systems where 
many different attacking rules and load patterns coexist, we need a 
better model, and a new level of description. 

From this perspective, applying real data of Beijing's traffic net- 
work, i.e. the subway network and the bus network, we find that the 
traffic system is robust for the NRb attacking model. Moreover, the 
efficiency of transferring people flow for removing a node of the bus 
network randomly (NRs) is higher with only increasing the tolerance 
parameter of the subway network as > ao than that of the bus 
network > ao. This means that the increment of some buses 
may release the traffic jam caused by NRs attacking model if the 
damage is limited. However, increasing the nodes capacity of the 
subway network as > ao will much sharper increase the efficiency 
to transfer people flow than that of the bus network. It is expected 
that the proposed model will reflect some important properties of 
real-life traffic systems when special properties of the subway are 
taken into account. 

Model. To illustrate how the flow redistributes on interrelated net- 
works, we introduce a network model. The network is considered as a 
weighted and undirected graph G with N nodes. The value of the 
weighted arc 6^ is a measure of the effect between adjacent nodes / 
and j. 

The load L/ on node / is the ''hetweenness' centrality, i.e., the 
probability of all the shortest paths passing through it^^'^^. The node 
capacity Ci by using the Motter-Lai model is deflned analogously as^^: 

Q = (l + a)L/, /=1,2, •••,iV. (1) 

The capacity parameter a e [0,1] is the control parameter represent- 
ing the tolerance capacity. If one node is removed from the depend- 
ency subway network (due to a failure or attack), it affects the 
shortest paths between some other nodes. Then the loads of many 
nodes are changed. If it is greater than its capacity, then it will over- 



load, and the efficiencies of its adjacent edges will decrease. Remote 
nodes will reselect its shorted path. It triggers further load adjustment 
over and over again until the loads of all remaining nodes are less 
than their capacities'^. 

With the above definitions, we can simulate the dynamic process 
of redistribution of people flows for the initial removal of a node in 
the related networks. It is well-known that all the sites or nodes on a 
line cannot work after one interdepended site is destroyed. Without 
loss of generality, we assume that a node / of the dependent network S 
is randomly selected as the initial damaged one. So the people at the 
same line of node / in network S will choose one of the vehicles in the 
connected network B (for more details, see Fig. 1). 

The redistribution rule produce a decrease of the efficiency (refer 
to Methods section) of the interrelated networks. In order to com- 
pare the effects of different attacks for the network's robustness 
against cascading failures, we adopt the following three basic attacks 
in the proposed cascading model: 

(1) Attack on the nodes randomly (NR): an easy attack strategy, 
here we compare the result of selecting a node of the dependent 
(NRs) and connected (NRb) network randomly and removing 
it; 

(2) Attack on the node with the highest load (HL): a common 
attack strategy, it is to remove the node with the highest load 
of the system. In heterogeneous networks, e.g., scale-free ones, 
HL is more likely to trigger cascading failures generally; 

(3) Attack on one node of the dependent chain with highest sum 
loads (HSL): it is to select the line number of the depended 
network in ascending order of the sum loads in the depended 
network and then to remove one node of the chain with the 
highest loads. 

Database in random network and real network. We analyze the 
robustness of different random networks with connectivity links and 
dependency chains. 

Network models. We study the robustness of ER and BA networks 
with connected lines and dependency chains. Here we consider two 
interrelated networks (cf. Fig. 2). Without loss of generality, the 
initial ER and BA networks are assumed to share the same size 
initially (i.e., the number of nodes N = 1000). Then choose edges 
from all possible ones of the ER network with the probability]?' = 
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Figure 2 | Interrelated random network model with connectivity lines 
and dependency chains. Suppose that there are six nodes in the initial 
network B at t = 0. Select two groups of nodes (Bi and B2 as chain 1; B4, B5 
and as chain 2) and add their corresponding nodes (Si, S2, S4, S5 and S^) 
at t = 1. Here the nodes of Si, S2, S4, S5 and constitute network S. 
Connect these new nodes of each group with the dependent chains as the 
new network S. Nodes Si and S2 are selected to be chain 1; nodes S4, S5 and 
are selected to be chain 2. Both networks are coupled by interconnected 
links of corresponding nodes from both networks, e.g., the arc of Bi and Si. 
tsy tBs and ts are the values of spending time passing two connecting nodes. 

0.01; the BA network grows according to the preferential attachment 
principle that each new node with m = 5 edges has larger possibility 
attaching to existing high degree nodes. Then select nodes by the 
following rule: for a Normal Gaussian distribution with an average 



number of nodes <s> = 10 in each dependent chains and the 
standard deviation cr = 0.5. The initial values of spending time 
passing two connecting nodes is ts = 2,ts= 1 and tss = 0.5 (eq.3). 

Real networks. Further we extend the analysis to real systems, the 
bus and subway networks of Beijing, and check the relationship 
between E and a, which are the average efficiency and the 
tolerance parameter. We choose all the 227 subway nodes and 924 
bus nodes of the lines which have 6 nodes in common with the 
subway nodes at least, and exert one of the three attack types on 
the interrelated networks. The largest load node of the systems is 
on the highest load line. For the NR attacking model, we give the 
average result over 10 trials. 

Results 

In Fig. 3 we find that the efficiencies of both the ER and BA network 
decrease with the increment of the percentages of dependency nodes 
(Pd). The efficiency of the ER network is about 0.03 when the per- 
centages of the dependency is 0 and the initial rate of moving nodes is 
0.5. However, for the same efficiency, the initial rate of removing 
nodes decreases to 0.3 if the percentage of the dependency increases 
to 0.6. Moreover, BA networks are a little more robust than ER net- 
works with the same increment of po- The efficiencies of BAs 
decrease slower than those of ER networks. 

In Fig. 4(a) we show the results of removing a node of the bus 
network randomly (NRb) for different values of the tolerance para- 
meter a in Eq.(l). We find that after both systems have reached a 
stationary state, the average efficiencies of the networks increase for 
increasing both as and a^. However, the average efficiency of the 
system with NRb for as = 0,1, = 0 is bigger than that for as = 
0, a^ = 0.1. The efficiency of the network is highest for as = 0.1, a^ = 
0.1. 
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Figure 3 | Efficiencies (eq.5) of ER and BA networks under randomly removing nodes (NR). p is the initial rate of removing nodes, p^ is the percentages 
of the dependent nodes, (a) Relations of the overload efficiency of the network, p and pj^ with nodes of the network removed randomly, (b) Relations of 
the overall efficiency of the BA network, p and p^ with nodes of the network removed randomly. The value of the tolerance parameter a is 0.03. 
Each of the networks has 1000 initial nodes. Each point is from an averaging over 10 realizations. 
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Figure 4 | Efficiency for the interrelated bus B and subway S network of Beijing, (a) and (b) are the processes of the redistribution of people flows in the 
interrelated networks with the attacking models of NRb and NRs respectively for different tolerance parameters ocs (the capacity parameter of the subway 
network) and (the capacity parameter of the bus network). The efficiency of the system for (star) increases (as =0) is lower than that as as 
(circle) increases (a^ =0) under both attacking models of NRb and NRs. The efficiency of the network with NRb attacking model is always greater than 
that with NRs attack strategy for the same tolerance parameters. For different effects of the dependent chains, there are huge error bars for NRs. 



We also study the average efficiencies of networks with attacking a 
node of the subway network randomly (NRs). In Fig. 4(b) the same 
situations were found. We also find the efficiency of the networks 
with NRb attacking model is always higher than that with NRs attack 
strategy. Huge error bars show there are diverse effects for the system 
attacking nodes of different dependent chains. 

In order to study the influence of the distributed flow on the 
efficiency of interrelated networks, we compare the efficiencies of 
these two networks with our four different attacking regimes in 
Fig. 5(a). We find that the efficiency of the system under the attacking 
model of NRb is highest. The efficiency under the attacking model of 
NRs is the next one. For the highest load node belonging to the 
highest sum loads chain of the subway network, the efficiency of 
the system under the attacking model of HL is almost equal to that 
under the attacking model of HSL for each value of a g [0,1]. Both of 
them are lower than that with NRb and NRs attacking strategies. 

We further study the efficiencies of networks under the envir- 
onment of attacking one node which randomly changes as increasing 
only the tolerance parameters of the subway network (as) or the bus 
network (a^) (see Fig. 5(b)). The efficiency of the system with the 
removal of a node of bus network (NRb) is higher than that of 
removing a node of the subway network for each value of ot e 
[0,1]. That is to say, removing a node of the subway network causes 
more serious damages than that of the bus network. The efficiency of 
the network as increasing ocs (ocb —0) is higher than that increasing ocb 
{ocs =0) for both random removing attacking models NRb and NRs. 
That means the efficiency of increasing the node capacity of the 
subway network will sharper increase the efficiency of the system 
than that of increasing the node capacity of the bus network. This 
confirms the results of Fig. 4. 

We further investigate the relationship among the tolerance para- 
meter a and the efficiencies of each network (see Fig. 6). In Fig. 6(a) 
we can see that the efficiency of the traffic system for NRb attacking 
model grows quickly to a constant for increasing the tolerance para- 
meter of the subway network (as) and the bus network (a^). 
However, it grows gradually for NRs (cf. Fig. 6(b)). The overall effi- 
ciency of the traffic system under the attacking model of NRb grows 



rapidly to recover the efficiencies with increasing a^ in a small region 
for some constant as and vice versa (Fig. 6(c)). So the traffic system is 
robust for responding to the attacking model of NRb. There exists a 
small critical threshold of the tolerance parameter ao for NRb. We 
find that there exist non-equilibrium phase transitions at the point of 
ao ~ 0.1. The efficiency of the network increases quickly with the 
increment of the tolerance parameter for a < ao, but it will recover 
immediately when a > ao. That is, small increasing of the tolerance 
parameter of the bus network could decrease the traffic pressure 
quickly. This explains why the pressure of the traffic congestion is 
relaxed by increasing some buses during rush hours. In Fig. 6(d), we 
find that there also exists a small critical threshold of the tolerance 
parameter ao ~ 0.1 for NRs. The efficiency of the network for as < ao 
and as = 0.1 is lower than that for ocb < ocq and as = 0.1. However, it 
is greater when as > ao and ocb > ocq respectively. Both efficiencies 
increase gradually with the increment of the tolerance parameter. So 
the efficiency of the traffic system recovers gradually as as and a^ 
increase. Increasing the number of buses can release the traffic jam if 
it is limited. But the traffic conjunction can be greatly slaked for 
increasing the tolerance parameter of the subway network (as) when 
a > ao (also could see Fig. 6(d)). 

Discussion 

Summing up, we have simulated a flow redistribution model of two 
types of interrelated networks, i.e. dependent network and connected 
network, and considered the special property of the subway's struc- 
ture that all the nodes in the same dependency chain are out of work 
once one node is failed. Within this framework, we have studied 
effects of the tolerance parameter of nodes and uncovered that the 
traffic system of Beijing is robust for an NRb attacking strategy. It will 
be best to transfer people flow if the tolerance parameter is larger 
than the critical thresholds. Moreover, we find that the efficiency has 
a non- equilibrium phase transition at the point of lower ao. This 
implies that increasing the number of small buses can substantially 
decrease traffic pressure during jam periods. Additionally, the effi- 
ciency will be recovered within a very small region of the tolerance 
parameter a. However, the damage of attacking the node of the high- 
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Figure 5 | (a) and (b) compare the efficiencies of the networks with four different attacking strategies, (a) is the efficiencies of the networks for four 
different attacking strategies as the tolerance parameters of both bus and subway networks increase at the same time. The green hues (circle) are for the 
traffic system under the attack of the highest load node (HL); the blue lines (square) are for the traffic system under the attack of the highest sum loads line 
of the subway network (HSL) ; the black lines (triangle) are for the traffic system when we remove a node of bus network randomly (NRb) , and the red lines 
(star) are for the traffic system when we remove a node of subway network randomly (NRs) . The efficiency of the system under the attacking model of NRs 
as the tolerance parameter increases are highest, following by that under the attacking model of NRs. And it is the same lowest with HL and HSL 
(b) is to compare the efficiencies with attacking models of NRb and NRs respectively for different tolerance parameters a. The red lines (square and star) 
are for the efficiencies of the networks as increasing the tolerance parameters e [0,1] when = 0. The black lines (circle and triangle) are for the 
efficiencies of the networks as increasing the tolerance parameters G [0,1] when as = 0. All data points of attacking model of Ni^ are averaged over ten 
times. 
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Figure 6 | In (a, b), the color 3D curved surfaces are the relations of the efficiency of the traffic system, changing the tolerance parameter of the subway 
network ( a 5) and the bus network ( a^) . ( a) is the relations of the overall efficiency of the traffic system, a 5 and with a node of the bus network removed 
randomly (NRb). (b) is the relations of the efficiency of the subway network subtracting the efficiency of the subway network removed randomly (NRs). 
(c) and (d) show some corresponding section planes of (a) and (b), respectively. The black lines (triangles) are the efficiencies of the networks as increasing 
the tolerance parameter of bus network (a^) when (Xs = 0.1. The red ones (squares) are the efficiencies of the networks as increasing the tolerance 
parameter of the subway network (ocs) when = 0.1. All the data points are averaged over ten times. 
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est load and the node in the highest sum loads line of the subway 
network is most serious for Beijing's traffic system. Furthermore, 
increasing the tolerance parameter of the bus network is more effec- 
tive than that of the subway network for a < ao- However, this is 
contrary to a < ao. These findings will be useful for relieving traffic 
congestion of our daily life. As the proposed model allows for effec- 
tive simulations, it could also be used for the analysis process of other 
interrelated networks. Moreover, BA and ER networks with depend- 
ency chains and connectivity links are studied. We find that increas- 
ing the percentages of dependency nodes could decrease the 
efficiency of the networks for the NR attacking strategy. The efficien- 
cies of a BA network is a little more robust than that of an ER network 
with the increment of the percentages of dependency nodes (p^). 

We emphasize the generic character of the cascading failures in 
Beijing's traffic networks; it can be analogously employed an other 
city traffic systems with dependency chains and connectivity links, 
such as Shanghai, Mumbai, Sao Paulo, etc.. Our model is readily 
applicable to these situations that all the nodes of the same depend- 
ency chain will be out of work if one node is attacked. The progress of 
cascading failures is not limited to assess the effect of overload nodes 
of city traffic systems, but is also relevant to many aspects of systems 
with dependency chains and connectivity links. However, this 
important aspect has received little attention so far. 

Methods 

First we show how the damage of a single node (refer to three attacking models in 
Model section) is sufficient to decrease the efficiency of the entire system simply based 
on the dynamics of redistribution of the flows on the networks. Suppose a node of 
network B is attacked at the beginning (see Fig. 1). This leads to changes of the most 
efficient paths (the shortest ones) between nodes and consequently the loads of nodes 
in network B redistribute, creating overloads on some nodes. Note that the overloaded 
nodes are not removed. Instead, for the efficiencies of the arcs, we adopt the following 
iterative rule at each redistribution time t> 2 (see Fig. 1): 



.(f+i)= 



:,(0) 



,(0)i 



iB,(0<CB, 



Otherwise 



(2) 



where it is supposed that the efficiency of the arc is equal to the initial value if the load 
of node / is not greater than the node capacity Cj. Otherwise, its efficiency will 
decrease meaning that the efficiency of all the arcs passing through node / are reduced 
so that eventually the flow will take alternative paths. 

We assume that the spending time passing through the connected nodes i and j, Up 

is inversely proportional to the efficiency of the arc, so it is defined by: tij ~ — .In this 

paper, we suppose the initial values of spending time passing a subway and bus 
network arc are tg = I and ts = 2, respectively. The initial value of time passing 
between a subway site and a near (less than half distance of two connecting nodes of 
the bus network) bus site is tss = 0.5. We also suppose that all users go along the 
shortest time path between two sites. Then we have 



1 



(3) 



where A;,, k2,... km are the nodes on the path from node / to j. The efficiency Eij of a path 
between all vertices i and j, can then be defined as the harmonic mean^° of the 
efficiency of the shortest time path between them 
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The robustness of the whole network is measured by the average transmission effi- 
ciency of the network: 



E(G)- 



1 



N(N-1) 



N(N- 



-1) 1 1 



(5) 



4. Zhou, D., Stanley, H. E., D'Agostino, G. & Scala, A. Assortativity decreases the 
robustness of interdependent networks. Phys. Rev. E 86, 066103 (2012). 

5. Parshani, R., Buldyrev, S. V. & Havlin, S. Interdependent networks: Reducing the 
coupling strength leads to a change from a first to second order percolation 
transition. Phys. Rev. Lett. 105, 048701 (2010). 

6. Huang, X., Gao, J., Buldyrev, S. V., Havlin, S. & Stanley, H. E. Robustness of 
interdependent networks under targeted attack. Phys. Rev. E 83, 065101(R) (2011). 

7. Lee, K.-M., Kim, J. Y., Cho, W.-K., Goh, K.-I. & Kim, I.-M. Correlated multiplexity 
and connectivity of multiplex random networks. New J. Phys. 14, 033027 (2012). 

8. Peng, H., Wei, N., Li, L., Xie, W. & Yang, Y. Models and synchronization of time- 
delayed complex dynamical networks with multi-links based on adaptive control. 
Phys. Lett. A 374, 2335-2339 (2010). 

9. Parshani, R. et al. Critical Effect of Dependency Groups on the Function of 
Networks. Proc. Nat. Acad. Sci. USA 108, 1007 (2011). 

10. Buldyrev, S. V., Parshani, R., Paul, G., Stanley, H. E. & Havlin, S. Catastrophic 
cascade of failures in interdependent networks. Nature 464, 08932 (2010). 

11. Parshani, R., Buldyrev, S. V. & Havlin, S. Interdependent Networks Reducing the 
Coupling Strength Leads to a Change from a First to Second Order Percolation 
Transition. Phys. Rev. Lett. 105, 048701 (2010). 

12. Brummitt, C. D., Lee, K.-M. & Goh, K.-I. Multiplexity- facilitated cascades in 
networks. Phys. Rev. E 85, 045102(R) (2012). 

13. Cardillo, A. et al. Emergence of network features from multiplexity. Sci. Rep. 3, 
1344 (2013). 

14. Tan, F., Xia, Y., Zhang, W. & Jin, X. Cascading failures of loads in interconnected 
networks under intentional attack. EPL 102, 28009 (2013). 

15. Bashan, A., Parshani, R. & Havlin, S. Percolation in networks composed of 
connectivity and dependency links. Phys. Rev. E 83, 051127 (2011). 

16. Youn, H., Gastner, M. T. & Jeong, H. Price of anarchy in transportation networks: 
efficiency and optimality control. Phys. Rev. Lett. 101, 128701 (2008). 

17. Li, G. et al. Towards design principles for optimal transport networks. Phys. Rev. 
Lett. 104, 018701 (2010). 

18. Gao, J., Buldyrev, S. V., Havlin, S. & Stanley, H. E. Robustness of a network formed 
by n interdependent networks with a one-to-one correspondence of dependent 
nodes. Phys. Rev. E 85, 066134 (2012). 

19. Gao, J., Buldyrev, S. V., Havlin, S. & Stanley, H. E. Robustness of a Network of 
Networks. Phys. Rev. Lett. 107, 195701 (2011). 

20. Gao, J., Buldyrev, S. V., Havlin, S. & Stanley, H. E. Networks formed from 
interdependent networks. Nat. Phys. 8, 40 (2012). 

21. Schneider, C. M., Araujo, N. A. M., Moreira, A. A., Havlin, S. & Herrmann, H. J. 
Mitigation of malicious attacks on networks. Proc. Nat. Acad. Sci. USA 108, 3838 
(2011). 

22. Schneider, C. M. et al. Suppressing Cascades of Load in Interdependent Networks. 
Proc. Nat. Acad. Sci. USA 109, E680 (2012). 

23. Liu, G. L., Peng, H. P., Li, L. X., Sun, F. & Yang, Y. X. Improving resource 
utilization in hierarchy network by optimizing topological structure. Eur. Phys. J. 
B 85, 63 (2012). 

24. Gomez, S. et al. Diffusion Dynamics on Multiplex Networks. Phys. Rev. Lett. 110, 
028701 (2013). 

25. Hu, Y., Ksherim, B., Cohen, R. & Havlin, S. Percolation in Interdependent and 
Interconnected Networks: Abrupt Change from Second- to First-order 
Transitions. Phys. Rev. E 84, 066116 (2011). 

26. Goh, K.-L, Kahng, B. & Kim, D. Universal Behavior of Load Distribution in Scale- 
Free Networks. Phys, Rev. Lett. 87, 278701 (2001). 

27. Latora, V. & Marchiori, M. Efficient Behavior of Small-World Networks. Phys, 
Rev. Lett. 87, 198701 (2001). 

28. Motter, A. E. & Lai, Y. C. Cascade-Based Attacks on Complex Networks. Phys. 
Rev. E 66, 0651 02 (R) (2002). 

29. Morris, R. G. & Barthelemy, M. Transport on coupled spatial networks. Phys. Rev. 
Lett. 109,128703 (2012). 

30. Crucitti, P., Latora, V. & Marchiori, M. Model for Cascading Failures in Complex 
Networks. Phys. Rev. E 69, 045104(R) (2004). 



Acknowledgments 

This work was partly supported by the National Natural Science Foundation of China 
(Grant Nos. 61100204, 61170269, 61121061), the Asia Foresight Program under NSFC 
Grant (Grant No. 61161140320), the Beijing Natural Science Foundation (Grant No. 
4142016) and the Beijing Higher Education Young Elite Teacher Project (Grant No. 
YETP0449), the Beijing Center for Mathematics and Information Interdisciplinary 
Sciences, and by the Government of the Russian Federation (Agreement No. 
14.Z50.31.0033). 

Author contributions 

Z.S. and H.P. designed research; Z.S., H.P., L.L. and J.K. performed research; Z.S., J.X., L.L. 
and Y.Y. explained the results; Z.S., J.K., J.X. and Y.Y. wrote and corrected the paper. 



1. Szabo, G. & Fath, G. Evolutionary games on graphs. Phys. Rep. 446, 97-216 
(2007). 

2. Castellano, C, Fortunato, S. & Loreto, V. Statistical physics of social dynamics. 
Rev. Mod. Phys. 81, 591-646 (2009). 

3. Barabasi, A.-L. The network takeover. Nat. Phys. 8, 14-16 (2012). 



Additional information 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Su, Z. et al. Robustness of Interrelated Traffic Networks to 
Cascading Failures. Sci. Rep. 4, 5413; DOI:10.1038/srep05413 (2014). 



SCIENTIFIC REPORTS | 4 : 541 3 | DOI: 1 0.1 038/srep0541 3 



6 



I This work is licensed under a Creative Commons Attribution-NonCommercial- 
ShareAlike 4.0 International License. The images or other third party material in this 
article are included in the article's Creative Commons license, unless indicated 
otherwise in the credit line; if the material is not included under the Creative 



Commons license, users will need to obtain permission from the license holder 
in order to reproduce the material. To view a copy of this license, visit http:// 
creativecommons.org/licenses/by-nc-sa/4.0/ 



SCIENTIFIC REPORT; | 4 : 541 3 | DOI: 1 0.1 038/srep0541 3 



7 



